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ABSTRACT 

• • 

The use of the random pKase approximation with exchange (RPAF) for 

cross 

calculating partial and total photo Ionization \ sections and photo- 
electron angular distributions for open-shell atoms Is examined for atomic 
chlorine. Whereas the RPAE corrections In argon (Z=18) are large, we 
find those In chlorine (Z=17) to be much smaller due to geometric factors. 
Hartree-Fock calculations with and without core relaxation are also 
presented. Sizable deviations from the dose-coupling results of Conneely 
are found. 


€ 
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I. INTRODUCTION 

If one excludes the lightest elements (i.e., 1<Z<10), then it is 
fair to say that the role of electron correlations in atomic photoioniza- 
tion processes has only been extensively studied for the rare gases and 
a few other closed-shell systems. The reasons for this are, firstly, 
the availability of detailed experimental rare gas photoabsorption spectra,^ 
and secondly, the relative theoretical ease of dealing with spherically 
symmetric targets having only a few final state channels. The key to 
the current theoretical understanding of these closed-shell photoabsorp- 
tion spectra has been the random phase approximation with exchange (RPAE). 

2 

The RPAE calculations of Aniusia, Cherepkov, and collaborators, and of 
3 4 5 

Wendin, Lin, and Starace have amply demonstrated the important influence 
that virtual excitations of pairs of valence electrons have on closed 
shell photoabsorption spectra and the necessity of taking these electron 
correlations into account in order to obtain good agreement with experiment. 
Consider specifically the history of some of the theoretical work concerning 
argon, one of the most thoroughly studied atoms. Hartree-Fock,®*^ intra- 

D 

channel, and close-coupling calculations failed to reproduce its experi- 
mental photoionization cross section Not until the RPAE calculation 
2 

of Ar»usia ^ al_. was good agreement (-10%) between theory and experiment 
obtained. Furthermore the importance of the electron correlations that 
the RPAE includes has been confirmed in argon by the very accurate many- 
body perturbation theory calculation of the cross section by Kelly and 
Simons. in addition, Burke and Taylor have recently done an R-Matrix 
calculation of the argon photoioni;idtion cross section that gives agreement 
with experiment comparable to that of the RPAE calculations. Thus while 
for very heavy closed-shell atoms it has been found necessary to take 
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into account core relaxation effects in addition to the electron corre- 

lations included in the RPAE, for the medium- heavy atom argon (Zsl8) 

the RPAE has given calculated cross sections that are in very good agree- 

11 1? 

ment with both many-body perturbation theory' and R-matrix' calculations 
and with experiment.^® 

Despite its success, the RPAE as originally developed for atomic 

14 2-5 

photoabsorption by Altick and Glassgold. and as used subsequently, 

has been largely limited to the theoretical treatment of closed-shell 

atoms. The only exceptions have been: (1) An extension of the '<PAE by 
15 

Cherepkov , to the treatment of photoionization by open-shell 

atoms having a half-filled valence shell (e.g., atomic nitrogen); and 

(2) An extension of the RPAE by Dalgaard^^ to the calculation of oscillator 

strengths for discrete transitions in open-shell atoms having either two 

electrons or two vacancies in the valence shell (e.g., silicon). Very 

recently, however, the general case has been solved and there are now 

two RPAE theories suitable for calculating the photoionization cross 

section of an arbitrary open-shell atom: these are tte theories of 
17 18 

Armstrong and of Rowe and Ngo-Trong . While these two theories differ 
in the requirements they place on the excitation operators that induce 
transitions between initial and final states, they are otherwise very 
similar. 

In this paper we use the open-shell RPAE theory of Armstrong^^ 

for the first time to explore electron correlation effects cn the photo- 

ionization cross sections and photoelectron angular distributions of 

atomic chlorine. We have chosen to study atomic chlorine (Z=17) since 

. < 

it is adjacent to argon (Z=18) in the periodic table and therefore might 
be expected to have similarly strong electron correlations that, as with 

- — ... — . . _ . _ V i d 
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argon, are accounted for adequately by the RPAE. In addition, it has 

been suggested that atomic chlorine, as well as the other halogen atoms, 

might in the near future be more readily studied experimentally than other 

open-shell atoms having Z>10. As further motivation for our study, we 

note again that at present the photoabsorption properties of approximately 

75X of the elements in the periodic table, mostly open-shell atoms, are 

largely unstudied either experimentally^ or theoretically (except in 

20 21 

the independent electron approximation). * 

Section II describes the reduction of Armstrong's^^ equations of 
motion to a f jrm suitable for numerical computation. The computation 
itself is described in detail. In brief, we have employed a basis set 
of discrete and continuum wavefunctions computed in the Hartree-Fock 
(HF) approximation and we have ignored coupling betw^-en photo?lectron 
channels. Section III presents our HF and RPAE calculations of total 
and partial photoionization cross secti^'^s and photoelectron angular 
distributions. We find that in contrast to argon, the electron corre- 
lations included in the RPAE are rather weak for atomic chlorine. We 
expect this result to hold even if future calculations include coupling 

between photoelectron channels. Comparison is made with a close-coupling 

22 23 

calculation of Conneely. * Section IV summarizes our results and 
conclusions concerning electron correlations in atomic chlorine and 
other halogen atoms. 

II. SOLUriON OF THE OPEN-SHELL RPAE EQUATIONS 

The equations-of-motion derivation^^’^^ of the RPAE starts by assuming 
a particular expansion for an excitation operator that induces transi- 
tions from the ground state (e.g., for photoionization this is the electric 
dipole operator). By taking appropriate commutators of this operator with 
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the Hamiltonian of the atom one obtains equations for 
the coefficients in the expansion of the excitation operator. These 

coefficients are all one needs to obtain numerical values for 

matrix elements of the excitation operator and hence for cross sections. 
17 

The Armstrong open-shell RPAE equations for the coefficient 
vectors X and Y of odd-parity, single-electron excitation operators are 
given most simply in matrix form: 



The matrix A represents the interaction betv/een singly-excited final- 
state configurations, while the matrix D represents the interaction between 
the ground state and doubly-excited states. The matrix U distinguishes 
Eq. (1) from the analogous equation for the closed-shell RPAE for which 
U becomes the unit matrix. U differs from the unit matrix because of 
the multiplicity of parent states .that an open-shell term level has. 
Finally, the excitation energy is labelled w, which equals the photon 
energy in a pliotoionization interaction. 

General expressions for the matrices A, D, and U in LS-coupling 
are given by ArmstrOivg^^ in terms of an arbitrary basis of independent 
particle wavefunctions. VJhile these general expressions appear to be 
quite complex, they simplify considerably upon making a few restrictions. 

In this paper we make the following restrictions: 

(i) We consider excitations from the open-subshell only. The 

27 

experience gained from using the RPAE for closed-shell atoms seems to 
indicate that inter-subshell interactions are important only for subshells 
having a small photoionization cross section relative to the cross section 
for nearby subshells. Neglecting interactions between subshells serves 
to make U diagonal and proportional to the square of a fractional parentage 
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thej 

coefficient for /ground state* open-subshell configuration. 

(ii) We choose for the excited orbitals an independent particle 
basis set of restricted HP wave functions calculated (in LS-coupling) 

in the field of the potential^® constructed from the set of HP dis- 
crete orbitals for the atom minus one of the open-subshell orbitals. 

This choice of basis set simplifies the matrix A by eliminating direct 

2 

intra-channel interactions betwaen continuum oroitals. 

(iii) We ignore any interaction between different final state 
channels, i.e., inter-channel coupling. To define these channels in 
LS-coupling, consider the photoionization of the lowest term level 
^^o^\ of the open-subshell nl^ of an atom,'^ by means of the electric 
dipole interaction: 








( 2 ) 


r 

We label the final state, on the* right in. (2), by ic where i, the channel 
index, is meant to specify particular values of S and L (trie spin and 

w V 

orbital angular momentum of the ionic core), of z (the photoelectron's 

orbital angular momentum’)', and of L (the total orbital angular momentum 

of the final state), c is the photoelectron energy, which may have 

discrete (i.e., negative) as well as continuum values. The electrostatic 

interaction is diagonal in L and hence neglect of inter-channel interactions 

amounts to neglecting matrix elements between states ic and je ' having 

differing photoelectron orbital momenta t. Close-coupling calculations 

indicate that these interactions do not have a large effect on the total 

29 

cross section except in resonance regions. We shall discuss our neglect 
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of the inter-channel interactions in more detail in Section III. Neglect 
of these interactions serves to make the matrix A diagonal in i and £ 
and equal to the product of the photon energy w and the diagonal matrix 

U. 

Expressions for matrix elements of U, A, and D appropriate to 
restrictions (i) and (ii) above are derived in Appendix A starting from 
the general expressions for these matrices given by Eqs. (27), (28), 
and (30) of Ref. 17. Restriction (iii) implies that we keep only those 
matrix elements, given in Eqs. (A. 10), (AJ.2), and (A.17)’ below, that 
are diagonal in the channel indices i and j. In order to indicate our 
method of solving the open-shell RPAE equations-of-motion in Eq. (1), 
it suffices to give here only the form of the matrices U, A, and D in 
the approximation that we make the restrictions (i)-(iii) above. From 
Appendix A, we have: 


^ie, jc' 


(3) 

^ic, je' 

4 (e+E^)U,, 

(4) 

r 

je' 

4 {..Dj ■ • 1 

10 le.je' 

(5) 


In Eq. (4), E^ is the binding energy in the ith channel. In what follows, 
then, we consider only interactions within channel i. the appropriate 
multichannel equations, however, may be easily inferred from the equations 
below provided one diagonalizes the matrix U.. at the outset of a multi- 

I J 

channel treatment. 

General Single-Channel Equations 


Solution of the single-channel open-shell RPAE equations-of-motion 
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starts by substituting Eqs. (3)-(5) in Eq. (1), noting that all matrix 


elements are real, and <(ividing both sides of to obtain 

coupled equations: 

the following 

^“"’''le. 1e' ' (“-=-El)*lcM 

(6) 


(7) 

where 


'^ie, ic' ~ ^ie, ic'^^ii* 

(8) 


and the integration over photoelectron energies e' includes a summation 
over discrete values of e'. Eqs. (6) and (7) can be uncoupled by first 
solving Eq. (7) for 




(9) 


and then substituting Eq. (9) in Eq. (6) to get a single equation for 




( 10 ) 


In Eq. (10) the matrix element H,. j.i(m) Is defined to be: 

I w I w 


^ic. ic'^***^ = ic'* 


( 11 ) 


Note that whereas the matrix element V. . , denotes the excitation of a 

I I c 

pair of electrons out of the open-subshell ground state to the final states 
ie and ie', the energy-dependent matrix element ic'(‘^) <*enotes an 


Indirect (i.e., second-6rder) Intra-channel interaction between final 
states i e and i c' . 

Solutions X. (od) of Eq. (10) may be written quite generally in the 
• e 

following form: 

' * (ur '^ic. 15<“> * H2) 

Here P indicates the Cauchy Principal Part is to be taken in any integra- 
tion over the singular denominator. K(co) is a reaction matrix, defined 
by the integral equation obtained upon substituting Eq. 12 in Eq. (10) 


i-cM = I'u.ieM * 




Note that X^g(w) depends only on the column of the reaction matrix K(u) 
having the index ic, where e is the photoelectron energy appropriate to 
the photon energy oj, 

e = o)-E^. (14) 

Lastly, B^(io) is a normalization factor determined from the open-shell 
RPAE normalization condition (Cf Eq. (33) of Ref. 17): 

/de{ = 6(u-u'). (15) 


In Appendix B we show that 
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Having obtained the coefficient vectors X. (w) and for an ar 

bitrary one-electron excitation operator we can now give the open-shell 

RPAE form for matrix elements cf the electric dipole operator in terms 

of the Coulomb interaction matrix elemtuts. Denoting the exact ground 

state by the bra (g| and the Hartree-Fock ground state by <g| and using 

similar notation for the final states ic, the reduced matrix elements of 

30 

the electric dipole operator are given by Ref. 17, Eq. (37) as: 


(g||r||ie)‘ 

(U,j)’’/de< 9 l|rl|ic> (X^^(ui) + Y,^(W) (17) 


Substituting Eq. (19) for Y..(w) in Eq. (17), then substituting Eq. (12) 

1 fe* 

for X^g(w), and making use of Eqs. (14) and (16), one obtains: 


(U)) 


(g||r||ic) * (<g||r||ie> * l’/de<g| |r| 


-/de:<g||r||ic >(^^) G,.^^j-(<a)) p + (vK-_-(o)))^)’‘‘, (18) 


where for compactness we have defined 


Gie.i-cM ^ ''ic.i-e ^ '’■''‘^’''iciCT^^^^T '^iC ,ic(“)- 


(19) 


Eq. (18) thus gives the open shell RPAE reduced electric dipole matrix 
element as the sum of three terms: (1) the zeroorder (i.e., HF) reduced 
dipole matrix element for photoionization from the ground state to the 
final state ic; (2) the contribution arising from indirect interactions 
between final states ic induced by virtual double excitations (cf. Eqs. 
(11) and (13)); (3) the contribution arising from the de-excitation of 


11 


the virtual, doubly-excited state (ic, ic) to the singly-excited state 
1 c by absorption of the .photon ou 

In sumnary, the reduced dipole matrix elements, from which cross 
sections and other measurable quantities may be calculated, are given 
In our one-channel, open-shell RPAE treatment by Eq. (18) in terms of 
zero-order reduced dipole matrix elements and in terms of two energy- 
dependent matrices K(uj) and G(w). i'(a)) is defined by Eq. (13) in terms 

of the indirect intra-channel interaction matrix l/(co), which in turn is 
defined by Eq. (11) in terms of the zero-order Coulomb matrix elements 

ic'* defined in terms of the matrix elements 

in Eq. (19). Any calculation of open-shell RPAE reduced dipole matrix 

elements thus starts from the computation of the matrix elements 

which are defined by Eq. (8) as the ratio of D. . , and U.., each of 

I e , 1 e II 

which are given in Appendix A. 

Calculational Details Specific to Photoionization 
of Atomic Chlorine 

Vie have considered the following photoionization reactions in atomic- 
chlorine: 

Ct3p®(^P) + hw Ct'*‘3p^(^V^L^) + e£(^L) (20) 

2S +1 

The allowed values of the ion term levels c L^, the photoelectron or- 
bital angular momentum t, and the final total orbital angular momentum 
L for the nine allowed channels i are given in the first three columns 
of Table I. Column four of Table I specifies the interaction matrix 

elements in each channel i in terms of the Slater integrals 

1 3 

R and R , where in atomic units 
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= R*^(3p3picJlE'£,) 


= /o/oI< P3p(«^l) Pci(^) Pap^^2) 

"< 

In Eq. (21) radial wavefunction of an n£ electron. 

The reduced dipole matrix elements are given by 


(21) 


<3p® ^P||r|13p^(^V\^)eA> 

“ a<3ptr!^^>^ 

• L 1 L Jl 

-(-1) '=3(2-5-(2wl).(2Ulf 

x(3p® ^P{|3p'’(^V\^). .p) <3p|r|£l> . 

where 

9 

<3p|r|e£> = Pe£(>") 


' 

I 

4 


{cl) 


(23) 


The coefficient a defined in Eq. (22) is given for each of the nine allowed 
channels in column five of Table I. The channel cross sections for 
photoionization of atomic chlorine art given in terms of the reduced 
dipole matrix elements .by 

0.^ = (0.4484 X 10"''V)(E.+e)<3p^ 1 |r ] |3p^(^V\^)eA>^ (24) 

31 

where the binding energies E. were taken from experiment and averaged- 
over fine structure levels; 
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j = 0.4768 a.u. (25a) 

E/, . = 0.5283 a.u. (25b) 

j = 0.6022 a.u. . (25c) 

Note that Eq. (24) requires E. and c to be in atomic units (1 a.u. = 

27,2108 eV) and the continuum radial wavefunctions Pg.j^(’“) to be 

normalized per unit energy (a.u.). 

The open-shell RPAE calculations employed a basis set of restricted 

HF excited orbitals computed in the V”” potential of atomic chlorine 

minus one 3p electron orbital. We shall call this basis the "unrelaxed 

HF basis. The potential was constructed using the discrete HF 

32 

orbitals for atcxnic chlorine given in the Clementi-Roetti tables. 

The first two discrete excited orbitals in each channel were canputed 

33 

using the multi -configuration HF computer code of Froese-Fisher. 

Continuum HF radial wavefunctions were computed by the method of Dalgarno, 

Henry, and Stewart^^using the computer code of S. T. Manson. These con- 

1 1 

tinuum HF orbitals correspond to those employed by Kelly and Simons 

and to the functions of Amusia et al 

An alternative HF calculation was also carried out using restricted 

N-1 

HF excited orbitals computed in the V potential constructed from dis- 
crete HF orbital s^^ appropriate to the alternative term levels of singly- 
ionized chlorine. This second HF basis, called here the "relaxed" HF 
basis, allows for core relaxation. Because of the slight non-orthogonality 
of the initial state (atomic) and final state (ionic) discrete HF or- 
bitals in this second calculation, we did not do an open-shell RPAE cal- 
culation using this basis set. The squared overlap integral between initial 
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and final state Slater determinant wavefunctlons is estimated from 

35 

similar calculations for argon photoionization to be better than 96%. 

We did not take this non-orthogonality into account and thus our 
“relaxed" HF results are too large by <4%. Our "relaxed" HF calculations 

7 ' 

correspond to similar calculations of Kennedy and Manson. 

In solving the open-shell RPAE equations all integrations over photo- 
electron energy e imply a summation over discrete energy values. This 
sunmation is iicmplished by first giving the discrete wavefunction for 

the last level nl included in our summation a continuum energy normal i- 

36 

zation by multiplying it by the factor 

n 

The summation over levels n't having n‘>n is then accomplished by starting 
the integration over e at rather than at threshold. 

In all, ve employed a basis set in each channel of 17 excited orbitals 
et. Two of these were discrete (either 3d and 4d or 4s and 5s) and 

t 

fifteen were in the continuum in the range 0£c<8 a.u. The various integral 

equations were reduced to a set of linear algebraic equations that were 

37 

solved by matrix inversion. 

Finally, note that all dipole matrix elements were calculated using 
the dipole length formula. Requiring the HF Hamiltonian to be gauge 

38 

invariant forces one to choose the length formula 'or HF calculations. 

2 

For RPAE calculations, Amusia show that length and velocity 

formulas are equal for closed shell atoms. Unfortunately we have not been 
able to prove this equality in the RPAE for open-shell atoms. 

In what follows we present zero-order results using both our "relaxed" 
and "unrelaxed" HF basis sets and open-shell RPAE results starting from 



the "unrelaxed" HF basis. 


111. RESULTS 
Partial Cross Sections 

In Tables II-IV we present the individual channel cross sections for 
photoionization of atomic chlorine. Each table presents cross sections 
for channels having in common a pa»ticular ion core term level. The last 
column of each table gives the sum of these cross sections, i.e., the 
partial cross section for photoionization of atomic chlorine leaving the 
ion in a particular term level. 

Comparing Tables II-IV with Table I we find the open-shell RPAE channel 
cross sections are smaller or larger than the "unrelaxed" HF channel cross 
sections depending on whether the coefficient of the Slater integral in 
the matrix element V- , is positive or negative. This is reasonable 

iE,ie ^ 

as seen by inspecting Eqs. (18) and (19), provided the term in Eq. (18) 
involving the matrix K(o)) has only a small influence on the reduced dipole 
matrix element. 

Since K- .■ ,(to) depends on the square of V. . , while G. . ,((o) 
depends linearly on V. . 1 * we can check the relative importance of 

1 C y 1 & 

these two matrices by investigating the magnitudes of the matrix elements 
'^ie ie'* maximum value of |V^ for any of the none channels and 

for any c or e' turns out to be 0.35 a.u. for the Ct ( S)ed ( D) channel. 

1 3 

This is due to the following circumstances: (a) the R and R Slater 
integrals are relatively large in this channel {e.g., they attain maximum 
values of 0.94 a.u. and 0,55 a.u. respectively); (b) the coefficients 
multiplying the R ' and R integrals (Cf. Table I) are both sizable and of 
the same sign in this channel. The other d-electron channels have peak 
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+ 1 2 

values of V- . . i that are 3G^i-90% of that in the Ct ( S)cd( 0) channel. 

1 C« I 1 w 

For these magnitudes of *Vj we find in all cases that the matrix 
elements of G(«i)) are at least an order of magnitude greater than those 
of K(u). Hence, it is the magnitude and sign of G. .• i(u) (which is 

I G 1 1 G 

nearly equal to V. . ,) that determines the observed behavior of the 

1 G j 1 c* 

open-shell RPAE channel cross sections relative to those of the "un- 
relaxed" HF channel cross sections. 

From Tables II-IV we see that in each channel the RPAE correction 
to the HF. cross section is frequently greater than 10%. The last column 
in each table, however, shows that summing the channel cross sections 
gives HF and RPAE partial ionic cross sections (i.e., corresponding to 
a particular ion term) that are usually well within 10% of each other. 

This reduction in the importance of the RPAE corrections for the partial 
ionic cross sections as compared to the individual channel cross sections 
is due to the alternating sign of the RPAE corrections in the individual 
channels which is discussed above. 

Figs. 1-3 present plots of the partial ionic cross sections in 
three different approximations. The solid lines are the RPAE cross sections, 
which Improve upon the "unrelaxed‘‘ HF cross sections, represented by the 
dashed lines. The dotted lines represent the "relaxed" HF cross sections. 

The relative behavior of the two types of HF cross section is similar to 
their behavior in irgon. In argon, the "unrelaxed" HF cross section®*^^ 

7 

is too high near the threshold, W'hile the "relaxed" HF cross section 
is more nearly the correct height, although shifted toward higher energies. 
The RPAE corrections for argon bring the "unrelaxed" HF cross section 
down toward the "relaxed" HF cross section near threshold. In chlorine, 
however, not only docs the RPAE cross secLiun show only small deviations 
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from the "unrelaxed" HF cross section, but in the case of the partial cross 
sections for the and ion terms the RPAE corrections to the ''un- 
relaxed" HF cross section widen the separation from the "relaxed" HF 
cross section near threshold (Cf. Figs. 2 and 3). If the "relaxed" HF 
cross sections near threshold are as reliable for atomic chlorine as they 
are for argon, then there would seem to be important relaxation effects 
in atomic chlorine that are not accounted for by the corrections in- 
cluded in the RPAE. Unfortunately the lack of experimental cross sections 
for atomic chlorine does not permit a judgment on whether the RPAE or 
the "relaxed" HF cross sections are more accurate near threshold. 

Photoelectron Angular Distributions 

We have calculated the photoelectron angular distribution asymmetry 
parameter 6 for photoionization leading to each of the three ion term 
levels. We use the angular momentum transfer formulation for $ in 
LS-coupling of Dill, Starace, and Manson.^^ jlie RPAE corrections to the 
"unrelaxed" HF phase shifts are given for each channel i by:^® 


’RPAE 


(w) = - arctan(7TK.£^.£((i))). 


(27) 


The maximum value of in any channel is 0.1 rad. The asymmetry 

parameters 3( P) and 3(*D) are plotted in Fig. 4; 3(*S) is plotted in Fig. 5. 
The cummulative effect of the differences between Hr and RPAE dipole matrix 
elements and phase shifts is no greater in tlie case of 6 than it is in the 
case of the partial ionic cross sections. Note that whereas Figs. 4 and 
5 show only the region of photoelectron 'energies 0<c<1.0 a.u., the near 
equality of the HF and RPAE results for 3 has been verified up to c=5.0 a.u. 


I 
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Total Photoionization Cross Section 

Fig. 6 presents our results for the total photo ionization cross 
section as a function of photon energy h\i. The difference between the 
"unrelaxed” HF and the RPAE calculations is even smaller here than it is 
for the partial ionic cross sections. This is because the RPAE correc- 
tions have alternate signs in the three ionic cross sections; for the 
P ion term (cf. Table II) the RPAE cross section is smaller than the 
HF one for hv»0.5 a.u.; for the ion term (cf. Table III) the RPAE 
cross section is larger that the HF one for 0,53 a.u.<hv<0.62 a.u,, 
but smaller for hv>0.62 a.u.; lastly, for the ion term (cf. Table IV) 

the RPAE cross section is larger than the HF one for h'^0.60 a.u. As 
is the case for the partial cross sections, here again substantial differ- 
ence exist between the "relaxed" and the "unrelaxed" HF cross sections 
that are not corrected by the RPAE. 

The crosses in Fig. 6 are the total close-coupling cross sections for 

22 23 

atomic chlorine of Conneely. ’ Our calculations do not include the 
inter-channel interactions that are included by Conneely. However, as 
discussed below, there is reason 'to believe that the large differences 
between our calculations and those of Conneely are not due to our neglect 
of inter-channel interactions. 

Comparison With Calculation of Conneely 

23 23 

The best comparison bet\»een our cross sections and those of Conneely ’ 
should be for the channel Ci-3p^(^P)+liv+Ct‘*’{ ^D)ed(^S). This d-electron 

C O 

channel can only interact with the s-electron channel Ct3p ( P)+hv^ 

Ct 3p^( S)es(^S), which should probably have only a very small effect on 
d-electron channel. Fig. 7 shows our results for this d-electron channel 


I 


I 
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and compares them with those of Conneely. In principle, except for the 
inclusion of the coupling to the s-electron channel, Conneely's cross 
sections for this d-electron channel (indicated by crosses) should be 
identical tc our "relaxed” HF cross section. For photoelec tjon energies 
e>0.3 a.u. our "relaxed" HF cross section is in fact identical to that 
of Conneely. For e<0.3 a.u., however, there are large discrepancies 
between the two calculations that most likely are not due to our neglect 
of coupling to the s-electron channel. 

It was easy for us to check the strength of the s- and d-electron 
channel coupling in argon and at the same time to verify that our "re- 
laxed" HF computational methods are working properly. We computed the 
"relaxed" HF cross section for photoionization of argon, 

Ar3p^(^S)+hv^Ar%^(^P)+e“(£=0, 2), 

# 

ignoring the coupling between the s- and d-electron channels. We compared 
our results to prior calculations of L. Lipsky^ which included (via close- 
coupling) the interaction between the s- and d-electron channels. We 
made sure to use exactly the same ^/*^”^ potential and the same binding 
energies as Lipsky and to take into account effects of non-orthogonality 
between initial and final basis sets. Thus the only difference between 
Lipsky's calculation and ours is the inter-channel coupling. Yet our 
results are within 4^^ of Lipsky's for all photoelectron energies in the 
range 0<e^.0 Ry. 

In summary, the discrepancy between our calculated cross sections 
for atomic chlorine and those of Conneely ni^ indicate that intcr-channel 
coupling is much more important in chlorine than in argon. However, if the 
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inter-channel coupling is in fact as weak in chlorine as it is between 
comparable channels in argon, then the discrepancy must have other, at 
present unknown, causes. Our agreement to within 4% with close-coupling 

9 

calculations of Lipsky for argon indicates that the discrepancies in 
chlorine are probably not due to any inefficiencies or inaccuracies in our 
computer code or method of procedure. 

Multichannel Open-Shell RPAE Calculation 

Because of the smallness of the intra-channel RPAE corrections to the 
single-channel cross sections, we have chosen not to proceed at this time 
to calculate the RPAE inter-channel interactions. Using the equations 
in Appendix A we have, however, calculated the angular factors for these 
interactions in order to obtain an estimate of their strength. The off- 
diagonal elements of the A-matrix give the direct interaction between 
photoelectron channels that is included in close-coupling calcjlations 
as well as in the RPAE. These off diagonal elements have a magnitude 
of =0.4 R^(3ped; e*d3p), which is not negligibl.i but at the same time not 
terribly large. The importance of these interactions can be determined 
precisely by comparison of the single-channel results we present in this 
paper with the results of close-coupling calculations using the same 
binding energies and potentials. Unfortunately the comparison between 
our results and the close-coupling calculations of Conneely is ambiguous 
at this point, as discussed above. 

Tlie off-diagonal elements of the D-matrix are only included in the 
RPAE and represent corrections to the inter-channel interactions arising 
from virtual double excitations. We find that the normalized off-diagonal 
elements Vj ^ i, defined by 

1 C , 1 E 
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''le.ic'S 




have angular coefficients of 0.1 for the R'(3p3pi ed e'd) Slater Integral. 
The magnitude of these off-diagonal angular coefficients is thus about 
one third the magnitude of those for the diagonal matrix elements (cf. 
Table I). Since the diagonal matrix elements have only a small effect on 
the single-channel cross sections, we thus do not expect the off-diagonal 
interactions. V. . , to alter the cross sections appreciably. 


IV CONCLUSIONS 


We have presented open-shell RPAE calculations of the partial and total 
photoionization cross sections and of the photoelectron angular distri- 
butions for atomic chlorine. Despite the fact that chlorine !tas only one 
f over electron than argon, the RPAE corrections in atomic chlorine have been 
found to be very much smaller than in ar,'*on. The direct reason for this 
is the smaller magnitude of the matrix elements that take into account 
the effect of virtual excitations of pairs of d-electrons from the ground 
state. In chlorine, in a typical d-electron channel, these matrix elements 
are =0.3 R^(3p3p;ed e’d), where is the largest Slater integral. In 
argon, by comparison, these matrix elements are =» 1.3 r\ Now in chlorine 
there ^ire six d-electron channels, whereas in argon there is only one. 

Thus one may argue that the RPAE corrections are, in fact, as strong in 
chlorine as in argon but are dissipated among six channels. Unfortunately 
the RPAE corrections do not have the same sign in all channels so that 
summing the individual channels to obtain the partial ionic cross sections 
or the total cross section reduces the difference between an RPAE 


calculation and a HF one. 
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It might be thought that perhaps only one of the eigenchannels of 
the photoelectron-ion interaction is strongly modified by the RPAE correc- 
tions, but that these modifications are reduced, in transforming to the 
experimentally observed channels given in Table I. To check this we 
diagonalized the A-matrix at a particular photon energy and transformed 
the D-matrix to this new representation. The transformed D-matrix had 
somewhat smaller off-diagonal elements and somewhat larger diagonal 
elements but unfortunately no single channel took unto itself all of 
the RPAE corrections represented in the D-matrix. The impossibility of 
simultaneously diagonalizing both the A and D matrices suggests then 
that the RPAE corrections are not concentrated in a single channel of some 
suitable representation. 

We are thus forced to conclude that the RPAE may not have as dramatic 
effects in open-shell atoms as in closed-shell atoms. Since this result 
in chlorine hinges on geometrical factor,- rather than dynamical ones 
(i.e., the Slater integrals are quite large in chlorine - only the angular 
factors are small), we expect similar results for the other halogen atoms. 
Alternative open-shell atoms for which the RPAE might introduce stronger 
correlations would probably have to be of one of the following kinds: 

1. Open-shell atoms having a reduced .number of photoionization 
channels, e.g., as for photoionization from ground or excited states 
having total orbital angular momentum L=0. Unlike the rare gases, however, 
these spherically symmetric atoms would still have several ionic term 
levels to photoionize to and thus there would still be more channels 

than for the rare gases. 

2. Open-shell atoms having RPAE corrections of the same sign in 
each photoionization channel so that the RPAE corrections produce a 
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cumulative effect on the partial and total cross sections. Alternatively, 
there might be some open-shell atoms for which the RPAE corrections are 
concentrated in only one channel. Unfortunately we know if no systematic 
way to pinpoint such atcmis, if indeed they exist at all. 

Lastly, it remains to be determined whether a proper kind of Hartree- 
Fock calculation (i.e., ''unrelaxed,'' "relaxed," or perhaps some other) 
is sufficient to correctly predict open-shell atom cross sections or 
whether there are other correlation effects, not included in the RPAE, 
which require a separate treatment. Unfortunately there is no experi- 
mental photoionization data as yet for atomic chlorine (or for other heavy, 

open-shell atoms) to answer this question even though the experiment for 

19 

atomic chlorine is certainly feasible. For this reason it is desirable 
to do other calculations for atomic chlorine, by different methods, par- 
ticularly because of the unresolved discrepancies between our results 
and those of Conneely. 
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Appendix A; Calculation of the U, A, and D Matrices 

We give here the expressions used to calculate the matrix elements 
A.j and D^jin the special case that we only consider excitations of a 
single ground state subshell nt^. We start in each case from the defining 
equations of Armstrong,^^ whose notation may be changed to that of this 
paper by making the following replacements: 

5. C. S'. C L^. S^.', L^' (A.la) 

S. L. S'.L* *^Sq. L^, Sq, L (A.lb) 

‘►et. c t . (A, Ic) 

The meaning of the symbols on the right of (A.l) may be inferred frcmi 

Eq. (2) of this paper. 

The general formulas of Ref. 17 for the U, A, and 0 matrices all 
involve the coefficients N and M (defined by Eqs. (17) and (18) of Ref. 17) 
which are in turn expressed in terms of multiconfiguration coefficients of 
fractional parentage (MCFP). Because they are ubiquitofi, it is convenient 
to evaluate the coefficients N and M first. 

The specific form of N which appears in this -particular calculation is 
N(ScLc* 9^q) * 

(A.2) 


where the CFP on the right are the HCFP defined by Armstrong, ' and g, y , 

c 

and Y indicate any additional quantum numbers necessary to specify the state 
?f the ground^ ionic core, and final configurations respectively. Use 
of Eq. (23) of Ref. 41 gives 


N 42 

where the CFP on the right is the. usual one for an nt" configuration, and 
T is the total number of electrons in the atom. The value of the second 
MCFP in Eq. (A. 2) depends on the specification of y* The most obvious 
and convenient form for y is 


(ySL) =1 (filled subshells) S L ,ctSL) (A.3b) 

in which case 


(ySI,{|YcS^Lc*‘='^> ' ■f’’* 


and thus 


N{SA,.Ye)l. 3U - ^ U? 9S L 


C 1 


(A.3c) 


(A.4) 


In what follows we shall not refer to the filled subiihells, as in Eq. (A. 3b), 
since they are passive. 

4 

The fora of M which appears in the calculation of the D matrix is 






1 (y3L{|y,'S^"L, 


M 



(A.5) 


*(Y'5,yiYc"Sc"Lc''> 


where the summation Is over y,". S,“, and L." and where [S ', L ...) 

C C V w W 
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Evaluation of Eq. (A, 5) requires specification of y* aod y^". We 
want y to have the form'givea in Eq. (A. 3b) and y' to be defined by 


.N-2 


WSoLq) = KC V2S2L2. (etc'OSjLjjS^L^). 


(A. 6a) 


The choice of y " is limited only by the requirement that it specify all 

N-2 

states of the configuration e£. Thp most obvious choice is 

With these choices of y, y' and the two MCfP in Eq. (A. 5) beccmie 
(■rSL{|Yc". Sj"Uc". «o) ■ (-1 ) V^c" V-c" 




", S , s "l^i* '■ ‘•cl , ^ ^ 

(‘0 4 4"j b 4 4" 


(A.7) 


and 


^ ^^0 ^^c^'c^l^o ^2^24’ ^0^ 


..V, I,-. L,,' P v| j. v| 


(A.8) 


where A , =2 if el“e't' and 4 i*! otherwise. Inserting Eqs. (A.7) and 

£ f C 

(A.8) into Eq. (A. 5) and summing over S " and L '' gives finally the desired 

t C 

expression for M: 

U/CU t B IlnlX- / 1 \L,.'+L,j+L +£ +S +-VS- '“So 
, yt^, yet) = (-1) c 3ooo‘*c 3 

*5 


c c 


xT-’((N-1)4^.^,)'‘ !S^'. L^-H S^. Lj. Sj]'^ 
J4' f4' ‘-3‘-c'li'4‘ 4 OPc' 4 4 ] 

(4 4 mV nh 4\h s h\ 
*(*■0 ' 444^14 7 e 44’ " o ^ 


(A.9) 
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U-Matrix 

• m 

For the special case of 'excitations from a partially filled shell 
to an empty shell. Ref. 17 shows that U is diagonal. The elements of U 
can be obtained either by combining Eqs. (30) and (32) of Ref. 17 or by 
substituting Eq. (A.4) above in the defination of U, given by Eq. (29) 
of Ref. 17. In either case one obtains 

^ic,jc' " ^^^c* **c* ^0* ^c * ^c * 


= 6(S^. S^') 6(L^, L^') 6(t. V) 6(e-e') N/f 

' 0^ 0 0 ' 0 'C C C 0 

^c 


(A.IO) 


A-Hatrix 

Using Eqs. (27) and (29) of Ref. 17 the matrix elements of A can be 
wri tten 

« 

^ic, je'~ ^b* ^*^c * ^c * ^ ^o* 




xN(S^’Lj‘. Y's'i'. 9*b). 


Substituting the form of N obtained in Eq. (A.4) gives 


{A.11) 


lE.je j2 i- 0 ^ 0 o ' 0 'c c c 0 


N-1 


Fc 


x(aI! g y's'i', ij 

'o’ oo'o 'c C C 0 


y,S^L^. et) SL|H|(t-'Y^'S^-L^'. c'f)SL) 


'C *1 * 


(A.12) 
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Evaluation of the A matrix is sin^lified considerably by use of a 
N»1 

y potential to calculate excited state wavcfunctions since in this case 

• (A.13) 

2S+1 

where E(YpS L , e£; SL) is the Hartree Fock energy of the- L state 

V w V 

N-1 

corresponding to an excited electron eS. moving in the V potential of the 
ionic level YpL^-S^* 

Since we have ignored interchannel interactions in this paper, the sum 

over Y* and y^' in Eq. (A.12) contains only one term: the one having the 
c ^ 

2S ■•■1 4 

particular term c of the ionic configuration 3p that is under 

consideration. Thus elements of A diagonal in the ionic core quantum 

numbers are given by: 

^c* ^0* ^c* ^c* ^0^ ^ (A.14) 

E(ScLc> SL) tf 

4 

Eq. (A.12) shows that matrix elements of A non-diagonal in the ionic 
core quantum numbers require evaluation of electrostatic matrix elements 
of the form 

i5j ttjKC\'Sc'Lc‘- SL). 

These matrix elements can be easily evaluated by the methods of Ref. 41 or 

43 44 

looked up in a number of published works. * 
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D-Hatrix 


The D-Matrix appropriate for this paper is given by Eq, (28) of 


Ref. 17; 




ie' “ *'c* ^0* ^c** **c'* %* 


> -H E •X'(S(.Lc’ ■'’*■ 0 ’ 

* N(Sc' L^', y't'V, gt|j) H(SjL^. v'tg, yeI)| 

* (^pWlI"19SoWl> 


(A.15) 


The matrix element of H in Eq. (A.15) is given by: 


41 


^YSo^olHlgVo^ = (N(N-1)/A^ (-1)^3 


Lo +1 


e.e 


^ ^V^o*’o^l^o T2S2L2» 


>- 3 ] 

(^0 ^0 ) 


(A.16) 


X u„l|c’'ll») (A„l|c''||i') ew'A') 


-QII^ I'-'' ‘'^0'' •' ' '~0 0’ 

. 41 45 

where the CfP is a two-body CFP., * The sum over y and y‘ in Eq. (A.15) 
is equivalent, in this case, to a sum over y^, y2^2^2' ^ 3 ^ 3 * 

stitution of Eqs. (A. 4), (A. 9), and (A.16) in Eq. (A.15) yields then the 
following expression for the D matrix elements: 

Di . , = (-l)^o‘^o'^'-c'^^c''^^''^V^ N(fM)/(2T^) 

ib,jb 

\S 3 L 3 f'-c- Sc* L 3 . S 3 . Lc* Sc)** (»oSo>-o{|"2'’^cScLc* *o> 


^2^24 


(A.17) 


,N-1 


^ ^^0 \^c4^l^0 ^24* ^o^^^oVo^l^o *^^244* ^o44^ 


N-2 







I 

I 


I 
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Appendix B; Normalization of RPAE Coefficients 

♦ • 

We wish to normalize the coefficients Y. (m) and X. (w) defined in 

1 w 1 W 

Eqs. (9) and (12) by finding the value of the coefficient B.(b>) in 
Eq. (12) that allows the RPAE normalization condition, Eq. (15), to 
be satisfied. We consider* each term in Eq. (15) separately. Using 
Eq. (12) the first term may be written; 

/deX^^(io) U^. =• Ujj Bt(io) B^(w') 


* - IsrsTT ^ 


+ p/dc K,^ ,-{■*>) ^ 


(w-e-E.) {a)'-e-E.) 




(Bl) 


Note that the photoelectron energies in Eq. (Bl) having a bar above 
them are fixed as follows: 

« 

e = «-E. and e' = w'-E. (B2) 

1 1 


Eq. (Bl) may be simplified by re-writing the product of .energy denominators 
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in the last term as follows: 


(«-£T E,)‘^ (w"-e-E^) 


-1 


(B3) 


) (a/-e'E^) (urc-E^) 


- — — } + ■/ 6(u)-u)') 6(a)- £ tE.) 

r \ * 


Substituting Eq. (B3) in Eq. (Bl) and rearranging terms gives: 
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/deX^^(u) • u:.Bt(u,)B.(uO 




X ^6(u-u') [1 


(B4) 


- Tm-c-rr 




il.ic-^) *t3^) '^u.ie-'<“'>j'’^ 


The two integrals in Eq. (B4) cancel identically. This may be 
verified by using Eq. (13) to substitute in the first integral 




and to substitute in the second integral 




(B6) 


L 
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Exact cancellation follows since V*{w) is Hermitian and real (cf. En.(ll)), 
Note that l/(o)) and both depend on u in Eqs.(B5) and (B6) even 
thnuyh the K'Hirttrix in Fq. (R5) depends directly on <.»' and that in Fq.(D6) 
on (I). This is because in considering the orthonormality properties of 
states with different energies we must require them to be solutions of 
the same Hamiltonian, i.e., the one containing V(u)), otherwise there is 
no reason for the states to be orthogonal if they have all other quantum 
numbers the same. 

Eq.(B4) thus simplifies to 

/deX+^{o.)U.,X, Ju') = iB^(u) I^U,., |^]«(.>-<a') (B?) 

A similar treatment of the second term in Eq.(15) leads to no such 
simplification. In particular, this second term is completely non-singular 
and no delta function 6(a> - (»'} can be factored out. Eq.(15) thus becomes 
upon substituting Eq‘. (B7): ■ 


6(o)-a)'') = iB^(u)) [1 + ^“(w) 


(B8) 


Taking the following integral of both sides. 


J O) + 

d(o' , 

oj - 6(0 (B9) 


results in 

which proves Eq.(16) In the text. 


(BIO) 
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TABLE I: Angular Momenta and Angular Momentum Factors for the 
Nine Allowed Photoionization Channels for Atomic Chlorine. 



1 

L 

'^le.le* 

a 


2 

2 

^ E.-. 
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2 

1 

CO 

o’* 

3p 

0 

1 

.1r1 
3 ^ 

+6** 


2 

2 

+1 r1 + 15_ r3 



2 

1 

1 pi 3 p3 

.s'* 

Id 

2 

0 

W ^ To 

-2/3** 

^0 

0 

2 


(10/3)‘ 

Is 

2 

2 

-L r’ - L 

T5^ 

- 2 / 3 ’“ 

Is 

0 

0 

2 r1 

(2/3 )'* 
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Tab’-e u; Cross Sections^ for Cil3p^(^P) + hv •*’ 
CH^3p*‘(^P)e)t( ^L) 


Photon Photo- 


Energy 

electron 

Energy 



ES 


^TOT 


hv(a.u. ) 

E(a.u. } 

HP 

RPA 

HP 

RPA 

HP 

RPA 

HP 

RPA 

O.U768 

0.000 

15.37 

ll*.69 

6.85 

7.59 

3.00 

3.10 

25.22 

rs.38 

0.5268 

0.050 

18.57 

17.53 

8.19 

9.07 

2.61 

2.69 

29.37 

29. 

0.5768 

0.100 

20. lO 

18.99 

8.62 

9.51 

2.27 

2.35 

31.29 

30.85 

0.6518 

0.175 

20.69 

18.93 

7.77 

8.53 

1.87 

1.93 

30.53 

29.39 

0.7268 

0.250 

19.09 

16.72 

5.92 

6.1*7 

1.55 

1.61 

26.56 

2l*.80 

0.8268 

0.350 

IU.36 

11.98 

3.1*5 

3.78 

1.23 

1.28 

19. oi* 

17.01* 

0.9768 

0.500 

6.95 

5.39 

1.25 

1.1*1 

0.893 

0.931 

9.09 

7.73 

1.2268 

0.750 

1.23 

O.O17 

0.161 

0.202 

0.569 

0.595 

1.96 

1.61 

1.1»768 

1.00 

0.089 

0.032 

0.003 

0.008 

0.395 

0.1*13 

0.1*87 

0.1*53. 

1.7268 

1.25 

O.Oll* 

0.031* 

0.016 

0.011 

0.291 

O.30I* 

0.321 

0.31*9 

1.9768 

1.50 

O.lll* 

0.139 

0.051 

0.01*5 

O.22I* 

0.233 

0.389 

0.1*17 

2.1*768 

2.00 

0.267 

0.270* 

0.098 

0.096 

0.ll*6 

0.151 

0.511 

0.517 

3.1*768 

3.00 

0.339 

0,301 

0.109 

0.113 

0.076 

0.078 

O.50I* 

0.1*92 

5.1*768 

5.00 

0.216 

C.201 

0.067 

0.071 

0.030 

0,030 

0.313 

0.302 


^Cross Sections given in units of 10 cni . 


^Cross Section for photoionization to the channel CJl '‘3p^(^)Ed(^D). 

^Cross Section for photoionization to the channel Cil*3p**(^)Ed(^P). 

d ^ ^ 3 _ 2 

”cross Section for photoionization to the channel CJl 3p ( r)£s( P). 

*Cross Section for photoionization to the ^ state of equal to the sum 

of the prior three columns. 
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Ttblc IV: Croso Sections** for C£-3p^(^P) hv 


Cl*3p**(^S)eA( ^L) 


Photon 

Energy 

Photo- 

electron 

Energy 



®T 0 T^ S) 

hv(a.u. ) 

c(a.u. ) 

HF 

RPA 

HF 

RPA 

HF 

FPA 

0.6022 

0.000 

7.80 

9.23 

0.436 

0.461 

8.24 

9.69 

0.6522 

0.050 

6.57 

7.28 

0.366 

0.387 

6.94 

7.67 • 

0.7022 

0.100 

U. 1«0 

1».70 

0.310 

0.328 

4.71 

5.03 

0.7772 

0.175 

2.13 

2.26 

0.245 

0.261 

2.38 

2.52 

0.8522 

0.250 

1.01 

1.10 

0.197 

0.210 

1.21 

1.31 

0.9522 

0.350 

0.378 

O.U36 

0.150 

0.161 

0.528 

0.597 

1.1022 

0.500 

0.078 

O.lOU 

0.104 

0.113 

0.182 

0 . 21 Y 

1.3522 

0.750 • 

0.000 

0.002 

0.063 

0.068 

0.063 

0.070 

1.6022 

1.00 

0.011 

0.007 

m 

o.o4i 

0.045 

0.052 

0.052 

1.8522 

1.25 1 

0.027 

0.023 

0.029 

0.032 

0.056 

0.055 

2.1022 

1.50 

0.039 

0.036 

0.022 

0.024 

0.061 

0.060 . 

2.6022 

2.00 

O.0I18 

0.048 

0.014 

0.015 

0.062 

0.063 

3.6022 

3.0c 

O.OU 

. 0.047 

0.007 

0.007 

0.051 

0.054 

5.6022 

5.00 

0.027 

0.029 

0.003 

0.003 

0.030 

0.032 


CL *l8 2 

“cross Sections given in xmits of 10 cm . 


^Cross Section for photoionization to the channel Cfc'*'3p (^S)efi(^D) 

/ 2 > 


*^Cross Section for photoionization to the channel C£,^3p**(‘^S)es(‘^S) 

*^Cross Section for photoionizaLion to the state of equal 
to the sum of the prior two columns. 
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FIGURE CAPTIONS 

1. Partial cross sections for photoionization of atomic chlorine 
leaving the ion in the state. Solid line: open-shell RPAE calcula- 
tion; dashed line: “Unrelaxed" HF calculation; dotted line; "relaxed" 

HF calculation. 

2. Partial cross sections for photoionization of atonic chlorine 
leaving the ion in the state. Solid line: open-shell RPAE calcula- 
tion; dashed line: "unrelaxed" HF calculation; dotted line: "relaxed" 

HF calculation. 

3. Partial cross sections for photoionization of atomic chlorine 
leaving the ion in the ^S state. Solid line: open-shell RPAE calcula- 
tion; dashed line: "unrelaxed" HF calculation; dotted line: "relaxed" 

HF calculation. 

3 

4. Photoelectron angular distribution asymmetry parameters B( P) and 
B(b). Solid line: open-shell RPAE calculation; dashed line: "unrelaxed" 

HF calculation. 

5. Photoelectron angular distribution asymmetry parameter 3(^S). Solid 
line: open-shell RPAE calculation; dashed line: "unrelaxed" HF calculation. 

6. Total photoionization cross section of atomic chlorine vs. photon 
energy. Ion term level threshold energies indicated by arrows. Solid 
line: open-shell RPAE calculation; dashed line: "unrelaxed" HF calcula- 
tion; dotted line: "relaxed" HF calculation. Crosses refer to the 
close-coupling calculation of Conneely. 

7. Partial cross section for photoionization of atomic chlorine leading 
to ejection of the photoelectron in the channel C2. 3p ('0)ed{^S). Solid 
line: open-shell RPAE calculation; dashed line: "unrelaxed" HF calculation; 
dotted line: "relaxed" HF calculation. Crosses indicate the close-coupling 
calculation of Conneely. 
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